The objective of this presentation is to give a brief overview of the theory behind the Damage Based 
Analysis (DBA) method including the derivation and practical application of the theory using the Python 
computer language. The Theory and Derivation will use both Acceleration and Pseudo Velocity 
methods to derive a series of equations for processing by Python. The results will compare both 
Acceleration and Pseudo Velocity methods and discuss implementation of the Python functions. Also, 
we will briefly discuss the efficiency of the methods and the amount of computer time required for the 
solution. In conclusion, DBA offers a powerful method to evaluate the amount of energy imparted into 
a system in the form of both Amplitude and Duration during flight and compare those results to the 
ground test qualification environments. Many forms of steady state and transient non-stationary 
vibratory motion can be characterized using this technique. DBA provides a more robust alternative to 
traditional methods such as Power Spectral Density (PSD) using a maximax approach. 
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Scope of the problem 


¢ Launch and Space Vehicles are subjected to complex, vibroacoustic non- 
stationary flight environments which vary in time during certain events in 
a mission profile. 

¢ Development, Qualification and Acceptance test specifications developed 
from historical data and analytical/empirical techniques don’t always take 
into account the non-stationary aspect of some portion of the flight data. 

¢ Flight data is assessed after each mission by comparing the flight results 
to prior qualification testing using a Maximax approach which can be 
potentially too conservative. 

¢ The need exists for a consistent, quantitative and accurate way to 
characterize and compare non-stationary flight data with stationary test 
environments. 


Historical Approach iS 


LAUNCH SERVICES PROGRAM = 


hs 


e Maximax Approach — Vibroacoustic Data 

— Typically, vibroacoustic flight data processed from instruments such as 
accelerometers are treated as a series of overlapping stationary time segments. 

— Maximax is a technique used to capture the Power Spectral Density (PSD) during 
each time segment. The goal is to obtain the maximum value of the spectrum at 
each frequency independently and loop through the total number of time 
segments. 

— The Maximax results are then compared to the test specification. 


e Limitations to Maximax Method 

— Peak responses in flight data can be considerably overestimated. 

— Maximax tends to result in excessive conservatism in the Qualification margin which 
may lead to more expensive designs and less reliability. 

— For arelatively immature design, component qualification may result in over-tests 
with excessive environments being applied which will impact cost and schedule. 

— Since Maximax is a purely mathematical transformation of the signal, results are 
dependent on the choice of processing parameters such as window size and 
overlap. 

— The physical properties of the system such as damping and fatigue from stress 
cycles are not considered. ; 
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Damage Potential 


Absolute Acceleration and Pseudo-Velocity Approaches 


Review of S.J.DiMaggio, B.H.Sako, S.Rubin , Absolute Acceleration and Scot I. 
McNeill Pseudo Velocity approaches. 


Absolute Acceleration Method, 
Damage Indicators, G1-G12 /% 
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Using Equations (4 & 5) in Appendix A, we can solve for the power spectral density (PSD) 
by plugging in «7 from eq(5) into eq(4) : 


7, oe 
= = (Sap ing) ae 


Eq(6) represents the power spectral density G,, which is a function of the max absolute 
acceleration response at each natural frequency represented by a SDOF system. Eq(6) is 
solved by computing an SRS at each frequency for a series of SDOF systems, recovering 
the Accel-Time responses and finding the maximum amplitude at each frequency. 


e Gl represents a Ty=60sec equivalent PSD input that ensures the test response or 
“SDOF base-drive response” of a component, envelopes the maximum amplitude 
response of a component in flight. 


e Note that G1 or amplitude damage indicator is related to maximum amplitude 
response of a component in flight and does not take into account fatigue. 


e Fatigue will be represented by damage indicators such as G4, G8 & G12 which are 
functions of the slope b for stress cycles on the S-N curve and are derived in 
Appendix A. Each of these indicators have a fairly complex mathematical form 
and will need to be derived each time the slope b changes, see Appendix A. 


e In order to illustrate the meaning of A2,,,/Q and G,, figures (1 & 2) are shown in 
Appendix A. 5 
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From Appendix B, It can be shown that the maximum stress o, is proportional to 
pseudo-velocity, o, = k o,, where o,, is the oscillator pseudo-velocity RMS. 


From eq(21) appendix A, we have: 


PVirax4 fn 
= Steal 21 
ot O in[faTol ee 


The resulting eq(21) represents the power spectral density (PSD), which is a function of 
the max pseudo-velocity amplitude response, PV,,,, at each natural frequency 
represented by a SDOF system. This result is similar to eq(6). 


From Appendix B eqs(23-25), the fatigue damage indicators G4 — G12 are: 


1 1 
a, = tthe [De G, = 18h *[ Oe iy = tthe *_ Oe 
SOM 8 Q 2Ap te Q | 720fnTo 


e Where D, is the number of fatigue damage cycles from flight data, counted using rainflow methods in a series 
of Amplitude Bins similar to the method used in the acceleration approach. 


Note that eqs(23-25) represent a much simpler form of the fatigue damage indicators 
compared to the acceleration method. Also, alternate fatigue damage indicators such as 
G6 or G10 could be easily derived for various values of slope b on the SN curve. 
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Practical Application - Python Implementation 


pyyeti.fdepsd.fdepsd( sig, sr, freq, Q, resp='absacce’, hpfilter=5.0, nbins=300, T0=60.0, rolloff=‘lanczos’, 


Fatigue Damage Equivalent — fdepsd python code ‘iS 
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ppc=12, parallel='auto'’, maxcpu=14, verbose=False) 


Parameters: 


Compute a fatigue damage equivalent PSD from a signal. See Appendix — C for more information 


sig (1d array_like) —- Base acceleration signal. 

sr (scalar) — Sample rate 

freq (array_like) — Frequency vector in Hz. This defines the single DOF (SDOF) systems to use. 
Q (scalar > 0.5) — Dynamic amplification factor @ = 1/(2Q)wnere Cis the fraction of critical 
damping. 

resp (sting; optional) — The type of response to base the damage calculations on: 


resp Damage is based on 
‘absacce’ absolute acceleration [1] 


‘pvelo’ —_ pseudo velocity [2] 


nbins (integer, optional) - The number of amplitude leveis at which to count cycles 
TO (scalar, optional) — Specifies test duration in seconds 


rolloff (string or function or None; optional) — indicate which method to use to account for the 
SRS roll off when the minimum ppc value is not met ‘linear 


maxcpu (integer or None; optional) — Specifies maximum number of CPUs to use. If None, itis 
internally set to 4/5 of available CPUs (as determined from 

multiprocessing -cpu_count Q) 

psd (2d ndaray; len(freq) x 5)-— The five columns are: [G1, G2, G4, G8, G12]: 


PSD Method Comparisons iS 
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Input Signal — based on arbitrary PSD input spec, converted to Acceleration-Time 


history. 


30-sec signal, sample rate=40kHz, 1.2M samples 


Input Signal - Specification Level 


Acceleration (g) 


15 20 
Time (sec) 


PSD Method Comparisons 


— PSD input spec 
— Welch PSD 

—  Maximax 

— Pseudo-Vel 


PSD (9°/Hz) 


Frequency (Hz) 


Acceleration (Gs) 


300 


Fatigue Damage Equivalent — mixed signal comparison iC 
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Fatigue Damage Equivalent (fdepsd) G1-G12 — Acceleration, Velocity methods 
compared to Maximax. 

Synthetic stationary Accel-Time signal injected with multiple non-stationary transient 
events approx. 200ms in length. fdepsd provides a better method to evaluate both 
combined steady-state and transient events in a composite signal. 


fdepsd test, Accel-Time History ‘93 fdepsd, ASD - Accel/Vel/PSD-Maximax methods comparison 
T T T T u T 


-200 


—300 


— PSD-G1-(Accel-method),Win64 et= 0.71 min 
- PSD-G1-(Vel-method),Win64 et= 0.53 min 

— PSD-G4-(Vel-method),Win64 

10*}) __  psp-G8-(Vel-method),Win64 

—  PSD-G12-(Vel-method),Win64 


PSD - Maximax 


ASD (G*2/Hz) 


1 1 


400 
0 


Time (sec) Frequency (Hz) 
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Performance considerations iS 
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For previous Fatigue Damage Equivalent processing which used Matlab like software, 
the average time to process the results for a signal > 300sec illustrated below was over 
45 minutes compared to 1.8 minutes for a Linux system. Note: Absolute Acceleration 
and Pseudo-Velocity methods for ASD computation correlate very closely. 


Total CPU time spent — processing fdepsd using Absolute Accel/Pseudo-Velocity 
methods on both Linux_64 and Windows-x64. Both systems set to 8-cpus, parallel 


processing. 


CPU time spent on Linux 64bit is very efficient and almost half that of Windows_x64. 
Most parallel CPU cycles on Linux consumed by user process with very little system/idle 


time. 


fdepsd test, ASD - Accel/Vel methods comparison 


ASD (G*2/Hz) 


\ | \ / 
Wy 
V 


PSD-G1-(AbsAccel-method),Linux et= 1.7862 min 
PSD-G1-(Pseudo-Vel-method),Linux et= 1.7944 min 
PSD-G1-(AbsAccel-method),Win64 et= 3.4526 min 
PSD-G1-(Pseudo-Vel-method).Win64 et= 3.3560 min 


Frequency (Hz) 


10? 
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Conclusions — 
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In summary, Damage Potential or Fatigue Damage Equivalent provides a powerful 
method for characterizing both the Stationary and Non-Stationary nature of a random 
vibration signal. Both the Amplitude and Fatigue characteristics of a signal are 
quantified based on characteristics of the system such as damping and stress cycles. 
Details for the damage indicators are summarized below: 


—  G,:A 60 sec equivalent PSD input that ensures the test response of a component 
envelopes the maximum amplitude response of a component in flight 


— G, Gg G,,: A 60 sec equivalent PSD input that ensures the fatigue damage of a 
component in test envelopes the corresponding fatigue damage in flight for a 
fatigue exponent of b=4,8,12. 


Pseudo-Velocity approach provides a more efficient method to apply fatigue damage 
exponents. 


The current Python implementation provides efficient Open Source software which is 
platform independent and can compute the results in a reasonable amount of time. 


Damage Potential Analysis provides a way to characterize the differences in energy 
imparted to a system during both Test and Flight and ensure that ground testing will 
envelope both the Amplitude and Fatigue qualification requirements of a component in 
Flight. 
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Appendix A — Absolute Acceleration Method 
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e Fatigue Damage Equivalent — FDE — Acceleration approach 


— The assumption can be made that the distribution of the peak values in a 
narrowband random vibration response follows a Rayleigh distribution for a given 
SDOF system, lightly damped stationary Gaussian input. 


— The Rayleigh probability density function for amplitude A, reference[1,5] is given 
by: 


AZ 
p(A) = (5) el-207) A>0 eq(1) 
e If we define the normalized and maximum amplitudes, see S.J.DiMaggio reference [1]: 


—2 
Amax Normalized max Amplitude response 
A*.ax Maximum Amplitude response 


e The maximum amplitude response, A2,,, is calculated from the Shock Response 
Spectrum (SRS) where the Acceleration, A2,,, is recovered from the Response 


Acceleration-Time history of the SRS. An ax = A* ../o* where o = (1-sigma) 
standard Deviation Accel(rms) response. 


SRS and Max. Response Acceleratigg 
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e SRS calculation 
— The SRS is based on the absolute acceleration transfer function response: 
wns 2 
_ a,(S)_ “9 *%n 
20) = ay(S) 524 One + coh eq(2) 
and the corresponding Frequency Response Function (FRF) Digital filter is: 
_ Bot B1z7*+ Boz~? 
HOY ata aa? eats) 


— The terms in equation(3) as defined in ISO 18431-4:2007(E) reference[4], represent a ramp 
invariant method or digital filter for solving the equations of motion for the SRS, refer to 
reference[4] for definition of FRF coefficients and more information. 


— Note: The coefficients in eq(3) are a function of : sampling freq. f, ,time interval T=1/f,, 
natural frequency f,,, natural angular frequency w,, and resonance gain Q. 


— From equations (22 & 5), S.J.DiMaggio reference[2] : 


Pirntise = Ara O° = 2InGfn To) eq(4) 
o2 = =G « f, *Q Miles Eq. eq(5) 


Where: o = (1-sigma) stand. Dev. Accel(rms) response, G = Accel. Spectral Density (PSD), 


fn = natural freq., T) = Ref. 60sec test time, Q = damping factor, a (zeta=damp. ratio) 
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Graphical representations, A2,,,/Q@ , G1 & G2 iS 
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Figure (1) reference[1] , represents a linear plot for 
the logarithm of time where response cycles exceed 
A*/Q. For T(A>) decreasing as a function of To, 


the function converges to 1/f, and the resulting 
2 


i AP 4 , . Az 
acceleration 9 is the maximum acceleration, —“* 


0.017) 


0.0017) 


Plot of G1 and G2, Freq = 664.3 Hz, Ty=382.4sec, A2.q,/Q=563.4 


1000 - 7 


Time (sec) - CycleCount/freq 


G1 


s Cycle Count (Time) vs (Amp*2/Q) 7 
G2 Figure (2) 
300 400 700 17 
Amplitude*2/Q 


Fatigue Damage Indicators: G4,G8,G12 
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In order to derive G4, we start with the definition of the Flight Damage Indicator D, as defined in S.J. 
DiMaggio reference [1], equation [14]: 


De = Q?/? S,CEPT(A,) eq(7) 
For b = 4: De = Q? Ehyray 
Cancel terms: De=* = ¥i,(A;)*T(A;) eq(8) 


In order to calculate the flight damage indicator in eq(8) for fatigue duration analysis, a rainflow 
process is used to count the cycles and place results into bins to account for the lower and upper 
bounds of the environment for each frequency, f,. Refer to reference[4], Standard Practices for 
Cycle Counting in Fatigue Analysis. 


Next, we need to derive the Test damage indicator D;. We start with the definition, from S.J. 
DiMaggio - reference[1], eq(18) : 


of eee a elect) a2) eq(9) 


=s5 


We can solve eq(9) by applying a series of substitutions and using multiple applications of integration 
by parts. The end result is: 


b=4_ 4 _ (~Amax ) —4 —2 
Dr To O [8 e ss (Artix +4* Amax + 8)] eq(10) 


Derivation of damage indicators — 
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If we equate D,; = D, for b=4 and solve for a2, we get : 
(=4maz) —4 —2 
yi(Ai)*T (Aj) = Tyo*[8 a ee (Amax + 4 * Amax + 8)] 
(A:)4 : 
o?= poe = Te8.Q eq(11) 


max.) _4 2 
2 (Amaxt+4*Amaxt8)] 


Solve for G = G4: 


Yi(Ad*T AD 
(~Amaz) _4 _2 
To [8 —e - (Amax+4*Amaxt8)] 


or G4 = 


G4 = APSD environment for fatigue damage indicator with exponent b=4. 
Ty = Time Slice — constant, default 60sec. 


_ 2 _ 
Amax = Ara O° = 2 In hi To) OF Amax = Amax 9 


Amax 18 the normalized max accel., (vector quantity) 


f, = Natural Frequency (Hz), evaluated for a series of SDOF systems spanning a specific 
Frequency range. 
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Derivation of damage indicators - Continued 
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ns = (Amax )? = (2 * tn me Toy* (vector quantity) 


T(A;) = Rainflow count of the Response Accel-Time History, recovered response from SRS of 
the input accel-time history, (vector) — Time Domain. See detailed counting of peaks and 
rainflow cycles in reference[4]. 


A, = Rainflow count of the Response Acceleration-Time History, SDOF base-drive response from 
SRS of the input accel-time history, (vector) — Time Domain. 


In a similar manner, 


(AN8T(CA; 
G8 = = 5 1/4 dy (AD eT (Aid) eq(13) 
p (=4max) _8 6 _4 _2 
T) (384=6. 2 /Aniey+8Amaet484 max t1924Amax 38h)! 
_ 2 1/6 yi {(Ai)12T (Ai) 
G12= ——* eq(14) 


(=4max) 12 —10 —8 —6 —4 22 
T> [384—e\ 2 (Amax+12Amax+120Amax+960Amax+5760Amax+23040Amax+46080)] 
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Appendix B, Pseudo-Velocity Method 
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Pseudo-Velocity Method i 
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The relationship between axial velocity and stress for a long thin rod iS dm, = PCVma,x where 
Vmax is the max. axial velocity. 
It can be shown that the maximum stress o, is proportional to pseudo-velocity, o, = k o,, where 
yy is the oscillator pseudo-velocity RMS. The corresponding pseudo-velocity transfer function is: 
G(s) = 2@& = ——*2 _ ; (pseudo-velocity)  eq(15) 


a(S) 52 4One cok 


and the corresponding Frequency Response Function (FRF) Digital filter is: 
H(z) = Bot Bixz t+ Bo*z~? 


1+ @1*Z~14+ Az*z~2 


; (pseudo-velocity) eq(16) 


For the Pseudo-Velocity approach, in a similar manner to the acceleration method, we define the 
Probability Density Function : 


P(PV) = aoe eq(17) 


where PV is the pseudo-velocity oscillator response. 


Integrating from arbitrary amplitude PV to infinity results in the probability that cycles have 
amplitude > PV or P(PV >). If the stationary environment duration is 7), the cumulative duration is 
T(PV >): 

-Pv2 


T(PV >) = Tp e2%bv eq(18) 
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Take the natural log of each side for eq(16), simplify terms where T(PV >) = > + as thé lower limit is 


PVAax 


reached for a where PV,2.,. is the maximum pseudo-velocity amplitude : 


In faTo] = eq(19) 


20p 


Recall from Scott I. McNeill paper, equation(15), reference[2], derived from Example 6.3, Force- 
Input/Disp.-Output system in reference[5] where the input is white noise : 


Ony = Paine where P, (f;,) is the PSD term, G,. eq(20) 
Substitute o,,, eq(20) into eq(19), rearrange terms and we have: 


= PVinax4T fn eq(21) 
Q In| fnTo] 


The resulting eq(21) represents the power spectral density (PSD), which is a function of the max 
pseudo-velocity response at each natural frequency represented by a Single Degree of Freedom 
(SDOF) system. This result is similar to eq(6). 


Calculate G4, G8 & G12, Fatigue Damage Spectrum — Frequency Domain, for the damage due to 
stress, eq(11), reference[2], we have : 


_s2 


D= “#T /*ns)stas where: p(S) = Se? eq(22) 


o% 


— where P(S) in eq(22) is the probability density function of the stress maxima, T is the total time 
of exposure to the stress environment, V,* is the number of positive cycles per unit time in the 
stress history, b is the fatigue exponent, c is a proportionality constant and S is the stress value ; 


of the peaks, see reference [2]. 
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Since maxima occur every (*/, ) seconds for a lightly damped oscillator response as illustrated in 
Figure(1), 4 = f,.. Wecan now rewrite equation (22) as : 
52 


— InT (© (1 +b) ve 
D= sc 40 S e295 dS 


Letu=S*% ; S= ju; du=2sds ; K= hf , alSo, for a property of the Gamma Function : 
Leds ae ; leta= ati perform integration, make substitutions and distribute sq.-root : 


D = #2 k(203,) "27 (; + 1) eq(22) 
Now substitute eq(20) into eq(22), Isolate o,,, and multiply by (2/b), cancel constants where 
c/ k® are (1), solve for P,(f,,) with b=4. Substitute G, = P, (f,) and r (5 + 1) =2,D=D,, 


T=T): 
Att fh Dr 
G4 = 
Q 2fnTo 


— Where D, is the number of fatigue damage cycles from flight data, counted using rainflow 
method in a series of Amplitude Bins similar to the method used in the acceleration approach. 


Similarly for Gg & Gy, where b = 8&12andIr (? + 1) = 24&720: 


Gy = fin *{ De eqi24) Gy = At | Pecans 
agp rer Q [720fTo 


eq(23) 
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Appendix — C, Python Code - Options 
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pyyeti.fdepsd.fdepsd(sig, sr, freg, Q, resp=absacce’, hpfilter=5.0, nbins=300, TO=60.0, rolloff=‘lanczos’, 
ppc=12, paralle/= auto’, maxcpu=14, verbose=False) 

Compute a fatigue damage equivalent PSD from a signal. 


e sig (7d array_like) — Base acceleration signal. 

e sr (scalar) — Sample rate. 

e freq (array_like) — Frequency vector in Hz. This defines the single DOF (SDOF) systems to use. 

e Q (scalar > 0.5) -— Dynamic amplification factor Q= 1/(2¢)where Cis the fraction of critical 
damping. 

e resp (string; optional) — The type of response to base the damage calculations on: 


resp Damage is based on 


‘absacce’ absolute acceleration [1] 


‘pvelo’ pseudo velocity [2] 
Parameters: 


e hpfilter (sca/ar or None; optional) — High pass filter frequency; if None, no filtering is done. If 
filtering is done, it is a 3rd order butterworth via scipy.signal.filtfilt(). 

e nbins (integer; optional) - The number of amplitude levels at which to count cycles 

e 10 (scalar; optional) — Specifies test duration in seconds 

e rolloff (string or function or None; optional) — Indicate which method to use to account for the SRS 
roll off when the minimum ppc value is not met. Either ‘fft’ or ‘lanczos’ seem the best. If a string, it 
must be one of these values: 


rolloff Notes 


‘fff ‘Use FFT to upsample data as needed. See scipy : signal : resample (). 26 


Fatigue Damage Equivalent — fdepsd python code “for inued 
d 


LAUNCH SERVICES PROGRAM s- 


ha 


‘lanczos’ Use Lanczos resampling to upsample as needed. See dsp. resample (). 


Apply a high freq. gain filter to account for the SRS roll-off. See srs.preroll1() for more 


Premer information. This option ignores ppc. 

dineae Use linear interpolation to increase the points per cycle (this is not recommended; method; it’s only 
here as a test case). 

‘none’ Don't do anything to enforce the minimum ppc. Note error bounds listed above. 

‘None Same as ‘none’. 


If a function, the call signature is: sig_new, sr_new = rollfunc(sig, sr, ppc, frq) . Here, sigis 
1d, len(time). The last three inputs are scalars. For example, the ‘fft’ function is (trimmed of documentation): 

ppc (scalar; optional) — Specifies the minimum points per cycle for SRS calculations. See also rolloff. Min. 
recommended ppc=10 which produces Max. errror at highest freq. of 8.14%, ppc=20, 2.05%. 

parallel (string; optional) — Controls the parallelization of the calculations: 


parallel Notes 


‘auto’ Routine determines whether or not to run parallel. 


‘no’ Do not use parallel processing. 
Use parallel processing. Beware, depending on the particular problem, using parallel processing can be 
‘yes’ slower than not using it. On Windows, be sure the £depsd() call is contained within: if __name__ 
== "— main __": 


e maxcpu (/nteger or None; optional) — Specifies maximum number of CPUs to use. If None, it is internally set to 
4/5 of available CPUs (as determined from multiprocessing. cpu_count () . 


e verbose (bool; optional) — |f True, routine will print some status information. 7 
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e maxcpu (integer or None; optional) — Specifies maximum number of CPUs to use. If None, it is internally set to 4/5 
of available CPUs (as determined from multiprocessing .cpu_count (). 

e verbose (boo/; optional) — \f True, routine will print some status information. 

e Arecord (SimpleNamespace class) with the members 

e freq (7d ndarray) — Same as input freq. 


Values Returned : 
e psd (2d ndarray; len(freq) x 5)- The five columns are: [G1, G2, G4, G8, G12]: 


Name Description 


The “G1” PSD (Mile’s or similar equivalent from SRS); uses the maximum cycle amplitude instead of 


ca the raw SRS peak for each frequency. G1 is not a damage-based PSD. 

G2 The “G2” PSD of reference [1]; G2 >= G1 by bounding lower amplitude counts down to 1/3 of the 
maximum cycle amplitude. G2 is not a damage-based PSD. 

G4, G8, ; ' 

G12 The damage-based PSDs with fatigue exponents of 4, 8, and 12 


e amp (2d ndarray; len(freq) x 5)- The five columns correspond to the columns in psd. They are the Mile’s 
equation (or similar if using resp='pvelo') SDOF (single DOF oscillator) peak response using the peak factor 
sqrt (2*log (f£*TO) ) . Note that the first column is, by definition, the maximum cycle amplitude for each SDOF 
from the rainflow count (G1 was calculated from this). Typically, this should be very close to the raw SRS peaks but 
a little lower since SRS just grabs peaks without consideration of the opposite peak. 

e binamps (2d ndarray; len(freq) x nbins) -— Each row (for a specific frequency SDOF 28 
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in binamps contains the lower amplitude boundary of each bin. Spacing of the bins is linear. The next column for 
this matrix would be amp[:, 0]. 


e count (2d ndarray; len(freq) x nbins) — Summary matrix of the rainflow cycle counts. Size corresponds 
with binamps and the count is cumulative; that is, the count in each entry includes cycles at the binamps amplitude 
and above. Therefore, first column has total cycles for the SDOF. 

e srs (1d ndarray; length = len (freq) ) — The raw SRS peaks version of the first column in amp. See amp. 

e var (1d ndarray; length = len (freq) ) — Vector of the SDOF response variances. 

e parallel (string) — Either ‘yes’ or ‘no’ depending on whether parallel processing was used or not. 

e ncpu (integer) — Specifies the number of CPUs used. 


resp (string) - Same as the input resp. 
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